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Stirling engine heaters are characterized by oscillating flows which affect heat transfer coefficients 
greatly. A 36-tube Stirling engine heater with a piston-link drive machine is used to study heat transfer 
characteristics of oscillating flows. The influences of the overall heating power, oscillating frequency and 
gas pressure on the heat transfer characteristics are investigated. To raise the working gas pressure is 
positive to decrease the wall temperature and to improve the heat transfer. When the pressure varies 
from 0.1 to 0.4 MPa, the wall temperature reduces 17 °C and the heat input increases 10 W. The time- 
averaged heat transfer coefficients reach the maximum value of 78.0 W/(m 2 K) among the testing condi¬ 
tions when the working gas pressure is 0.4 MPa and the revolution is 420 rpm. An oscillating heat transfer 
correlation is derived based on the experimental data which are obtained under situations close to an 
actual Stirling engine’s flow region, where Re and Re^, are in the ranges of 740-4110 and 12-71, respec¬ 
tively. The estimated uncertainty for the heat transfer coefficient is usually within 7.69%, and the deter¬ 
mination coefficient of a regression fitted correlation is 0.97. The proposed correlation is respected to 
predict the heat transfer coefficients of oscillating flows for practical design of tubular heaters, while 
the classical unidirectional steady correlations are not suitable, especially at high Re conditions. 

© 2013 Elsevier Ltd. All rights reserved. 


1. Introduction 

Stirling engine heaters are characterized by oscillating flows 
which affect the heat transfer coefficients greatly. A lot of studies 
on the flow and heat transfer characteristics of oscillating flows 
have been done to get a deep understanding of performance of 
Stirling engines. 

Simon and Seume [1] summarized the features of oscillating 
flows and heat transfer in Stirling engines, and found out that 
many modern Stirling engines were operated in or near the lami¬ 
nar to turbulent transition region. Table 1 lists some classical lam¬ 
inar and turbulent steady unidirectional correlations which are 
used to get approximations of convective heat transfer coefficient 
in Stirling engines. Urieli and Berchowitz [2] introduced a so-called 
Quasi steady flow model, as the Eq. (d), assuming that there was a 
steady flow condition in each integration time. Kuosa et al. [3] 
compared four laminar oscillating correlations, i.e. Eq. (e-h), and 
indicated that some of the laminar oscillating correlations gave 
almost the same Nusselt number to turbulent unidirectional equa¬ 
tions. However, many experimental researches showed different 
results. Shin and Nishio [4] adopted the Womersley number as a 
dimensionless value of oscillating frequency to consider the vis¬ 
cous effects, and derived a heat transfer correlation, as the Eq. (i), 
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of oscillating flows in tubes using a numerical study under an 
assumption of isothermal and extremely thin walls. Fan et al. [5] 
employed this correlation to predict heat transfer coefficient of a 
long tube of 2.1 mm inside diameter and 0.41 mm wall thickness. 
His experiments showed the correlation gave reasonable results 
at low oscillating frequencies and provided under-predicted values 
at higher frequencies. Kanzaka and Iwabuchi [6,7] focused on the 
influence of piston phase difference on the heat transfer coeffi¬ 
cients in a singular tube, proposed a correlation considering the 
piston phase difference, as the Eq. (j), which was verified later by 
the experiments using a Stirling engine. Bouvier et al. [8] experi¬ 
mentally investigated oscillating flows in a circular tube of 5.7 m 
long, 48.3 mm outer diameter and 3.2 mm thickness, and found 
that temperature oscillations in the wall could be ignored. Akdaga 
and Ozguc [9] studied the oscillating flow of water in a vertical 
annular liquid column with a constant heat flux, and obtained an 
empirical correlation, i.e. the Eq. (k), for a cycle averaged Nusselt 
number as a function of kinetic Reynolds number and dimension¬ 
less amplitude. These studies indicated that oscillating flow factors 
have great impacts on heat transfer performance, including kinetic 
Reynolds numbers, frequencies and fluid displacements. While 
there is no widely recognized conclusion, especially in correlations 
for heat transfer of oscillating flows. Besides, most of the above 
experimental studies were carried out using singularly slender 
tubes, not real Stirling engine heaters which always had tens of 
tubes. 
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Nomenclature 


A heat transfer area, m 

A 0 , A w , A r dimensionless fluid displacement 
c p working gas specific heat capacity at constant pressure, 

J/(kgK) 

c f correction factor for gas 

C correction term 

d tube diameter, m 

/ Darcy resistance coefficient 

f r Reynolds friction factor 

g gravitational acceleration, N/kg 

Gr the Grashof number 

h heat transfer coefficient, W/(m 2 K) 

H characteristic length, m 

I current, A 

L overall length of a tube, m 

n number of the heater tubes 

Nr revolutions per minute of the flywheel, rpm 
Nu the Nusselt number 

P pressure, MPa 

Pr the Prandtl number 

R 2 determination coefficient of the regression fitted corre¬ 

lation 

Re the Reynolds number 

Re max the maximum Reynolds number 

Rea, the kinetic Reynolds number 

T temperature, I< 

u working gas velocity based on mean piston velocity, m/s 

U voltage, V 

w uncertainty 

W 0 the Womersley number 

x max fluid displacement in heater tube, m 


Greek symbols 

a coefficient of cubical expansion, 1 /K 

s emissivity of material 

A conductivity of material, W/(m K) 

A dimensionless fluid displacement 

/u viscosity of fluid, N s/m 2 

v kinetic viscosity of fluid, m 2 /s 

cp heat, W 

cd oscillating frequency, rad/s 

Subscripts 
b buffer 

c cooler 

C convection 

cal calculated results by formerly proposed correlations 

e effective 

exp experimental results 

g working gas 

h heater 

i inner wall 

ins insulations 

o outer wall 

obj stands for insulations or heating medium 

oil heating medium 

p displace piston 

r regenerator 

R radiation 

sur the surrounding 

w tube wall 


This paper investigated the heat transfer characteristics of oscil¬ 
lating flows of a 36-tube Stirling engine heater, whose operating 
conditions were close to the flow regions of real engines, and 
derived an empirical oscillating heat transfer correlation to predict 
oscillating flow heat transfer performance for promoting practical 
design of tubular heaters. 

2. Experimental setup and methods 

An experimental apparatus was constructed to study the heat 
transfer characteristics of oscillating flows in a tubular heater. A 
schematic of the apparatus is shown in Fig. 1. The apparatus con¬ 
sists of a device body, an oil-bath heating block to heat the tubular 
heater, a heat sink via cooling water, a gas supplement unit, and a 
piston-link drive machine with a crosshead. The device body has 
three heat exchangers, i.e. a heater, a regenerator and a cooler, 
and the structure is modeled after displacement type Stirling 
engines. A motor connected to a crankshaft drives the displacer, 
which makes working gas flow cyclically between the heater and 
the cooler through the regenerator. 

The heater is made of an array of 36 U-shaped stainless steel 
tubes. These tubes form a sort of cage around and above the cylinder 
of a 110 mm diameter. The outer diameter and overall length of 
each tube are 6 mm and 298 mm, respectively. The annular wire 
mesh regenerator is set up in a stainless steel holder. The cooler con¬ 
tains 30 aluminium tubes with fins, whose inside diameter and 
length are 5 mm and 70 mm, respectively. The volumes of the hea¬ 
ter, regenerator, cooler and displacer sweeping are listed in Fable 2. 

At the start of each experiment, the seals were checked. The 
sealing units on the displacer contain three piston rings and a 


guide ring made of Rulon, preventing gas short-circulating the heat 
exchangers. Stepseal in combination with Variseal were chosen as 
piston rod seals to seal-off the pressurized working gas within the 
cylinder. O-rings were applied at other matching parts as static 
seals. Pressures in the cold end of cylinder and in the buffer were 
measured to test sealing effects of rod seals. A required initial pres¬ 
sure (0-0.8 MPa) in cylinder was controlled by the gas supplement 
unit. 

Nitrogen was applied as the working gas, and the revolution 
varied between 150 and 500 rpm during all experiments. The tem¬ 
perature of dimethicone was between 70 and 250 °C in the heat 
transfer tests. The cooling water rate was kept around 40 L/h. Test¬ 
ing data was recorded by a computerized data acquisition system, 
with measuring points and locations illustrated in Table 3. 

The mass of working gas fluctuates as gas travels through differ¬ 
ent temperature sections in a cycle. The reductions of the mass of 
working gas in Fig. 2 show the leakage behaviors at different rev¬ 
olutions. It is interesting to observe that higher oscillating fre¬ 
quency helps to suppress gas leakage. This is mainly because the 
seals expand due to friction generates more heat at higher 
frequency. Fig. 3 shows the static leakage behaviors at different 
buffer pressures. It is found that the smaller pressure difference be¬ 
tween the cylinder and the buffer, the better sealing. 

Dimethicone was used as the heating medium, i.e. an oil bath, 
which was heated by an electric ribbon heater and surrounded 
by thermal insulations. In heat transfer tests, the overall heating 
powers were maintained at a steady state and were evaluated by 
the voltage and current values. The temperatures of the outer wall 
of insulation, the heating medium, and the surroundings, were 
measured by seven thermocouples. The effective heat transferred 



Table 1 

Heat transfer correlations in literatures. 
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into the working gas was determined by the overall heating power 
minus heat dissipations, which is given by 10], 


(Pe UI ((Pins,C (Pins.R ) Poil.C ^ PoiLR ) 


( 1 ) 


(Pobj.C — h-objAobj (T obj Tsur) 


( 2 ) 


(Pobj,R = 5 . 67 s obj A obj [(To/jj/lOO ) 4 - (T sur / 100) 4 


(3) 


h-obj — {NUobj S-obj) / Uobj 


(4) 


Ntt-obj — 0.5 9 (Gr objPf obj) 


1/4 


(5) 


GV 0 bj — §^objU 0 bj(Tobj T sur )/V obj ’ 


1.43 x 10 4 < Gr <3 xlO 9 . 


( 6 ) 


The heater tubes immersed in the heating medium were instru¬ 
mented with thermocouples near the ends. The temperature distri¬ 
bution was assumed fairly uniform. Besides, the hypothesis of 
neglecting temperature oscillating in the wall is acceptable in this 
certain operating conditions [8 . Thus the average temperature of 
outer tube wall was detected and the average temperature of inner 
wall could be calculated according to the physical thermal conduc¬ 
tivity properties of the heater. Thermocouples were inserted into 
connections between the hot end of cylinder and the heater and be¬ 
tween the heater and the regenerator, respectively, to get the aver¬ 
age gas temperatures in the heater tubes. The time-averaged heat 
transfer coefficient was given by, 


T Wfi = T 


W,0 


(Pe(do d-i) 

2Xnd 0 L 


(7) 


hg — (p e / n diT(T W j T g ) 


( 8 ) 


3. Results and discussions 

3.1. Heat transfer experiments of The heater 

Fig. 4 details the heat transfer performance at various overall 
heating powers. Apparently, a higher overall heating power yields 
higher wall temperatures of the heater tubes. The effective heat in¬ 
put from the heat block to the working gas increases with the over¬ 
all heating power at first and then drops to some extent at the end 
when other operating conditions keep constant. Fig. 5 shows that 
the temperatures of the wall and the heating medium decrease 
when the working gas pressure raises and the overall heating 
power remains unchanged. Fig. 5 also indicates that raising the 
working gas pressure is positive to improve the effective heat in¬ 
put. When the pressure varies from 0.1 to 0.4 MPa, the wall tem¬ 
perature reduces 17 °C and the heat input increases 10 W. 

Figs. 6 and 7 present relations between the heat transfer coeffi¬ 
cient, the oscillating frequency, and the gas pressure. The time- 
averaged heat transfer coefficient is between 39.9-78.0 W/(m 2 K) 
during the tests. A higher frequency brings the oscillating flow 
with a larger core and a thinner boundary layer, which enhance 
the heat transfer. Fan et al. [5] reported that these characteristics 
could enhance the heat transfer both between the inner tube wall 
and the working gas and along the axial direction of the fluid. An 
improvement on heat transfer by adding the pressure in cylinder 
indicates that large mass flow rate is also beneficial. 

The accuracy of the volt-ammeter was taken as 1% ±2 divisions, 
and that of the K-type thermocouples could be 0.1 °C within the 
experiment range. The uncertainty analysis of the time-averaged 
heat transfer rate was performed with the main source of errors 
[9,15] as follows. 
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1-heater tubes, 2-cylinder, 3-displacer, 4-regenerator, 5-piston ring, 6-cooler, 7-frame plate, 8-piston rod, 9-piston rod seals, 
10-crosshead, 11-connecting rod, 12-crankshaft, 13-buffer, 14- ribbon heater, 15-thermal insulations 
I- device body; II- oil bath device; Ill-cooling water system; IV- gas supplement unit; V-drive mechanism 

Fig. 1. Schematic of experimental apparatus. 


Table 2 

Volume of each component. 


Item 

Volume/cm 3 

Heater 

70.76 

Regenerator holder 

202.32 

Regenerator matrix 

30.78 

Cooler 

42.23 

Displacer sweep volume 

134.30 


Table 3 

Measuring items and location of measuring points. 


Object 

Location 

Symbol 

Measuring 

instruments 

Pressure 

Cold end of cylinder 

Pc 

Pressure 

transducer 


Buffer 

P b 

Pressure gauge 

Temperature 

Hot end of cylinder 

T h 

K-type 


Interface between heater and 

Ph,r 

thermocouple 


regenerator 

Interface between cooler and 

Tc.r 



regenerator 

Cold end of cylinder 

Pc 



Outer wall of heater tubes 

Tw,o 



Immersed in dimethicone 

Poil 



Attached to the outer layer of 
insulation 

Pins.o 



Atmosphere 

Psur 


Revolution 

Flywheel 

Nr 

Tachometer 

Overall heating 

Ribbon heater 

U,I 

Ampere-voltage 

power 



meter 



Fig. 2. Leakage behaviors at different revolutions. 


Wh g — yj(dh g /d(p e ) + (dh g /dT W j) + ( dh g /dT g ) Wj s (9) 

The deviation of the calculated effective heat input is less than 
12.5 W. The estimated uncertainty for the heat transfer coefficient 
is 5.42-7.69%. Experimental results are listed in Table 4. 

For further investigations, test data was converted into dimen¬ 
sionless parameters. The gas flow inside the test section with a 
constant overall volume is driven by the displacer. Assumptions 
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Fig. 3. Static leakage behaviors at different buffer pressures. 


250 


200 
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CD 
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150 ^ 
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CD 

100 ^ 


50 


50 100 150 200 250 300 350 400 

Overall heating power(W) 


Fig. 4. Effective heat input and average wall temperature of the heater tubes vs. the 
overall heating power. 






Thus the Reynolds number based on the cross-sectional mean 
velocity and the kinetic Reynolds number for the flow in the heater 


could be obtained as: 


Re = udi/v 

(11) 

Re m = cod 2 /v 

(12) 


During the tests, the Re is in the range of 740-4110, and the corre¬ 
sponding Rea, is 12-71. Taking oscillating effects into consideration, 
the heat transfer correlation in the heater could be expressed as fol¬ 
lows. The correlation based on experimental results is regression 
fitted with the determination coefficient of R 2 = 0.97, seen in Fig. 8. 

Nu = 0.0162Re° ) 4 X 85 (13) 


Fig. 5. Effective heat input and average temperatures at different gas pressures. where, Aq - X max /dj is a dimensionless fluid displacement. 


are made that the flow field is one-dimensional and the velocity in 
cylinder is in accordance with the displacer movement. The work¬ 
ing gas cross-sectional mean velocity in the heater is written by 

[ 6 ], 

u = (d p /di) 2 X (2 x max Nr/n x 60) (10) 


3.2. Comparisons between experimental and theoretical results 

Ratios of the experimental Nusselt number and calculated value 
by some theoretical correlations in Table 1 are plotted in Fig. 9. If 
considering the deviation of ±20% is acceptable, the calculated re¬ 
sults from the Eqs. (c) and (i) are reasonable. It is noticed the Eq. (c) 
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Table 4 

Cycle-averaged experimental data and uncertainty analysis. 


Items 

co (rad/s) 

Rea, 

Ao 

U/(W) 

<Pe (W) 

Wcp e (W) 

Tyv,in (°C) 

T g (°C) 

M W/(m 2 K)) 

Uncertainty (%) 

Nu 

1 

21.36 

36.18 

181.05 

130.9 

100.96 

in 

100.45 

79.82 

55.00 

7.63 

5.54 

2 

21.68 

33.91 

181.05 

159.6 

118.97 

8.4 

121.87 

96.24 

52.15 

7.06 

5.07 

3 

22.51 

26.69 

181.05 

337.5 

230.63 

12.5 

205.70 

160.31 

57.10 

5.42 

4.89 

4 

33.51 

56.04 

181.05 

130.9 

100.14 

7.7 

99.14 

82.41 

67.26 

7.69 

6.73 

5 

21.68 

33.91 

181.05 

161.5 

126.63 

8.5 

121.87 

96.24 

55.51 

6.72 

5.39 

6 

29.85 

47.86 

181.05 

161.5 

125.14 

8.5 

112.81 

91.05 

64.64 

6.80 

6.35 

7 

33.51 

53.18 

181.05 

161.5 

124.28 

8.5 

114.44 

93.23 

65.84 

6.84 

6.44 

8 

39.79 

62.68 

181.05 

161.5 

123.51 

8.5 

114.70 

94.78 

69.67 

6.89 

6.79 

9 

43.98 

71.28 

181.05 

161.5 

128.60 

8.5 

107.37 

88.85 

78.01 

6.61 

7.70 

10 

33.51 

12.93 

181.05 

161.5 

114.04 

8.5 

131.76 

99.62 

39.87 

7.46 

3.85 

11 

33.51 

25.99 

181.05 

161.5 

116.58 

8.5 

125.81 

98.38 

47.77 

7.29 

4.62 

12 

33.51 

39.28 

181.05 

161.5 

118.22 

8.5 

120.55 

96.61 

55.48 

7.19 

5.39 

13 

33.51 

53.18 

181.05 

161.5 

124.00 

8.5 

114.44 

93.23 

65.69 

6.86 

6.42 



Fig. 8. Experimental Nusselt number as a function of Re w and A 0 . 



Fig. 9. Ratios of experimental and calculated Nu vs Re. 


gives Nu values generally smaller than the experimental ones. The 
Eq. (i) over-predicts Nu values a little at low Res, while under-pre- 
dicts at high Res. The other three equations, unidirectional turbu¬ 
lent correlations, give much bigger Nu for most Res, and 
deviations become larger when Re increases, which was also no¬ 
ticed by Kuosa’s [3 . 

These experiments confirm that the heat transfer characteristics 
of the oscillating flows are very different from those of the 
unidirectional steady flows, especially when Re is high. The kinetic 


Reynolds number, the Womersley number and the dimensionless 
amplitude are important parameters for the oscillating flows. 

4. Conclusions 

A 36-tube Stirling engine heater was used to investigate the 
heat transfer of oscillating flow characteristics. The structure of 
the experimental apparatus and the testing conditions are similar 
to the operating conditions of real Stirling engines. The experimen¬ 
tal results show that to increase the working gas pressure is a po¬ 
sitive way to improve the effective heat input and to decrease the 
wall temperature of heater tubes, reducing the requirement on 
materials. The time-averaged heat transfer coefficients between 
the inner tube wall and the working gas increase with the working 
gas mass flow rate which is relative to the pressure and the oscil¬ 
lating frequency. An empirical oscillating heat transfer correlation 
is derived based on the kinetic Reynolds number and the dimen¬ 
sionless amplitude, i.e. Nu = 0.0162Re° J 4 °Ao 85 , where Re and Re^ 
are in the ranges of 740-4110 and 12-71, respectively. The esti¬ 
mated uncertainty for the heat transfer coefficients is usually with¬ 
in 7.69%, and the determination coefficient of a regression fitted 
correlation is 0.97. The experimental results are also compared 
with formerly proposed correlations, indicating that the classical 
unidirectional steady correlations are not suitable for oscillating 
flows, especially at high Re conditions. The empirical oscillating 
heat transfer correlation provides time-averaged heat transfer 
coefficients of oscillating flows, expecting to be valuable for practi¬ 
cal design of tubular heaters. 
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